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The number of bound-state solutions of the Schrodinger equation for the screened Coulomb potential 
(Yukawa potential), -(Or) exp(-ar), occurs frequently in theoretical discussions concerning, for example, gas 
discharges, nuclear physics, and semiconductor physics. The number of bound states is a function of (C/a). 
Three upper limits for the number of bound states associated with the Yukawa potential are evaluated and com- 
pared. These three limits are those given by Bargmann, Schwinger, and Lieb. In addition, the Sobolev inequality 
states that whenever (C/a) < 1.65 no bound state occurs. This agrees to within a few percent of the numerical 
calculations of Bonch-Bruevich and Glasko. The Bargmann and Lieb limits and the Sobolev inequality are 
substantially easier to evaluate than the Schwinger limit. Among the three limits, the Schwinger limit gives the 
most restrictive limit for the existence of only one bound state and, therefore, is the best one to use for the ap- 
proach to no binding, i.e., 1.65 < (C/a) < 1.98. The Lieb limit is the best among the three when (C/a) > 1.98. 
The Bargmann limit is the least restrictive. 
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1. Introduction 

Many theoretical topics in such areas as plasma, nuclear, solid state, and semiconductor physics [1-6] 1 re- 
quire knowledge about the number of bound-state solutions of the Schrodinger equation for two particles 
that interact according to the Yukawa or screened Coulomb potential. For example, as semiconductor 
devices become smaller, required dopant ion densities and their spatial gradients increase. The density of 
states for such highly doped materials is needed to predict device performance and degradation mecha- 
nisms. The density of states involves summation over bound states and integration over continuum states. 
The disappearance of bound-impurity levels in semiconductors due to free-carrier screening of the Coulomb 
field of impurity ions is thereby fundamental in understanding the behavior of electronic devices. [7] 

During the course of research on the effect of donor impurities on the continuum electronic states in semi- 
conductors, the author has used the Bargmann limit [8], the Schwinger limit [9], and the Lieb limit [10] to 
evaluate upper limits for the number of bound states associated with the screened Coulomb potential. The 
analytic methods used here to calculate these limits are more general and may be applied to other fields 
such as those mentioned above. 

We shall show that 1) when the Yukawa potential permits only one bound state to exist, Schwinger's 
extension of Bargmann's limit gives a more restrictive limit for the number of bound states than do the 
Bargmann and Lieb limits; 2) when many bound states exist, the Schwinger and Bargmann limits agree to 
within a few percent and are not as restrictive as the Lieb limit; and 3) the Sobolev inequality [10,11] 
predicts no binding when (C/a) < 1.65. The limit 1.65 from the Sobolev inequality agrees to within a few 
percent of the limit 1.68 from the numerical work of Bonch-Bruevich and Glasko. [5] 
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2. Yukawa Potential and Schrodinger Equation 

The Schrodinger equation for two particles that interact according to the Yukawa potential is 

[V 2 +^-{£-- 7(r)}]iKr) = f (1) 

where E is the energy of the two-particle system, m is the reduced or effective mass, r is the relative coordi- 
nate vector between the two particles, and r = |r| is its magnitude. The spherically symmetric potential 
energy V(r) is assumed to have the Yukawa form: 

V(r) = -Cr'expi-ar). (2) 

For many applications, the Yukawa potential is an idealized one which gives a reasonable description of the 
physical phenomena under study. Expressions for the parameters C and a in terms of physical quantities 
then depend upon the application. For donors such as phosphorus or arsenic in silicon, the quantities C and 
a become, respectively, C = (e 2 /e) and a. = \lr,. The charge of the electron is e, the dielectric constant is e, 
and the screening length is r t . The last depends in part upon the temperature and the dopant density. It 
gives the extent to which free carriers screen the Coulomb field of the donor ion. 

The radial wave equation obtained from the Schrodinger equation, eq. (1), for the spherically symmetric 
potential V(r) is 

~§T Mr) + [^§-{E- V{r)} - *' + *> ] Ur) = , (3) 

where 

lMr) = Ea tm f (r)r |m (M). 

The quantities, E^, are constants. The normalization and orthogonality conditions for the spherical har- 
monics Yt„ (0,0) are 

f d<f> \ x dd sin Y hmt (0,0) Y hm2 (0,0) = b hh 6 mim2 . (4) 

In the following sections, we shall calculate upper limits for the total number of bound-state solutions iVof 
the Schrodinger equation, equation (1), and the number of bound-state solutions n t o( the radial wave equa- 
tion, eq. (3), for a given L Counting multiplicity, we perform the following summation to obtain the total 
number of bound states for a given spherically symmetric potential V. 

N[V] = f f o (2?+\)n ( [V), (5) 

where ( max is the largest value of £ for which n ( ^ 0. 

3. Bargmann Limits 

When the integral, 

h[V]= \; r \V(r)\dr, (6) 

is finite, the Bargmann formulation [8] gives the inequality that 

(2f + l)n, < 1 B [V\. (7) 
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Thus, for the Yukawa potential, eq. (2), we have 

I B [V]=(C/a). (8) 

From the inequality (7), we note that the maximum £ = t max for which n t =£ in the sum of eq. (5) has the 
upper limit 1/2 (I B — 1). Thus, we have the following inequality [9] for spherically symmetric potentials; 
namely, 

N<±-I B (I B + 1). (9) 

4. Schwinger Limits 

Schwinger [9] extends the Bargmann derivation to treat angular and spin-independent potentials and to 
give the number of states that lie at or below some chosen energy. His result for the upper bound I s to the 
total number of bound states N associated with a general spin-independent three-dimensional potential 



where 



N < Is , (10) 



c^v Wiiyi . 



(4tt) 2 ' ' |r-r 



We shall now proceed to evaluate the double integral in eq (1 1) for the case of the Yukawa potential given 
by eq. (2). Because the Yukawa potential is spherically symmetric, we expand the denominator in terms of 
the spherical harmonics; namely, 

T^T" = 4 * S i- e£i> 7r r " iW r '~ m • (,2) 

where r, and r^, are, respectively, the lesser and the greater values of |r| and |r'| and where 

The Legendre polynomials P™ satisfy the relation 

P ,- W =<-,)- ffzE| PrW . (14) 

Substituting relation (12) into eq. (11) and using for the volume element in spherical coordinates d 3 r = r 2 
sin dddd<t>dr, we perform the integrations over 6, </>, 6', and 0'by frequent reference to eqs. (4), (13), and (14) 
and obtain the result that for spherically symmetric potentials V,(r) 

h[V,\ = J V dr j V 2 dr' | Vs(r)\ \ Vs{r')\ 

(15) 



x f 1 r§L 

' =0 (2f+l) rf 2 



Equation (15) becomes for the Yukawa potential (2) 

/,(Y) = (C/af f ^L_ [ \; 2 exp(-z) dz 



(16) 
X( £ V 2Ul exp(-2r,) d V + £ v - 21 - 1 exp(-z7j) dr,\] , 
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where we made the substitutions z = ar, y =ar\ and then rj = (y/z) to obtain equation (16) from equation 
(15). 

Let us consider the first term in the curly brackets in eq. (16). Interchanging the order of integration, we 
obtain 

rfi(0 = lo V 2U1 (Jo "* exp(-z) exp(-z7/) dz} dr, 

We now examine the second term in the curly brackets of equation (16), namely, 

m = lo"* exp <-*) (Jfi,-^ 1 exp(-«,) d V ) dz . (18) 

Interchanging the order of integration in equation (18), we write the second term 

m = Ev^a+if)- 2 ^. (i9) 

By making the substitution rj = (1/x), it follows that d 2 (i) = d^l). 

From the above, we obtain the Schwinger upper limit for the total number of bound states associated with 
the Yukawa potential 

h[V] = (C/a) 2 f 2(2f+l)- 1 dffl . (20) 

Interchanging the order of integration and summation in eq. (20), 

Is[V) = (C/af Jo(l-77)- 2 L 2r ] 2e+l (2e+iy i 



and observing that when 77 < 1, 



/„ (1+m - r 2q 2fM 



we write eq. (20) in the form 

I S [V] = (C/af £ ln\(l+v)l(l- V )) dr] (21) 

We show in Appendix A that the integral in eq. (21) equals 1/2. Hence, our final result is 

Is[V] =\{Ib[V]V. (22) 

5. Lieb Limits and Sobolev Inequality 

In three dimensions, Simon [10] gives a bound for large (C/a) due to Lieb of the form 

iV< fl3 | \V(r)\ 3 ' 2 d 3 r = J L [V], (23) 

where a 3 < 0.116. For the potential (2), the Lieb limit h[V] becomes 

I L [V] = 0.116 x 6.062 [h[V]} 3n . (24) 

When (C/a) is small enough, bound states do not occur. The Sobolev inequality [10,11] predicts that 
bound states do not occur whenever hi[V] < 1, or equivalently, whenever 
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1 > 0.078 X 6.062 [I B [V]} 3 ' 2 m I sl [V] . 
That is, according to the Sobolev inequality, N = whenever (C/a) < 1.65 

6. Conclusions 

We have shown in sections 3, 4, and 5 that the upper limit for the total number TV of bound states 
associated with the Yukawa potential is for the Bargmann limit, eq. (9), 



for the Schwinger limit, 



and for the Lieb limit 



N<jI,[V]Ub[V] + D=n b , 



N<y [I B [V]) 2 = N s , 



N < 0.7032 {Ib[V]} 3 ' 2 = N L 



Because < N S [V] < N B [V] for all I B [V] > 0, we conclude that the Schwinger upper limit, N S [V], is a 
more restrictive limit than N B [V]. However, for large enough values of I B [V], the fractional difference [(N B 
— N S )IN B ], between the Schwinger and Bargmann upper bounds approaches zero as [\II B [V ]) approaches 
zero. 

When (C/a) > 1.98, N L < N s and the Lieb limit is best for large (C/a) and for many bound states. When 
(C/a) < 1.98, N s < N L and the Schwinger limit is best for the existence of only one bound state and for the 
approach to no binding, i.e., 1.98 > (C/a) > 1.65, and finally, the Sobolev inequality gives N = whenever 
(C/a) < 1.65. This value of 1.65 agrees to within 2 percent of the value (C/a) < 1.68 which Bonch-Bruevich 
and Glasko [5] determined numerically to give no bound state. 

And finally, we observe from table 1 that even the Lieb limit is more than a factor of 3 greater than the 
numerically determined number of bound states given in reference 5 for values of (C/a) > 8.92.* 

TABLE 1. Comparison of the Number of Bound States N Determined 

Numerically by Bonch-Bruevich and Glasko' and the Three Upper 

Limits of Bargmann [Eq. (9)], Schwinger [Eq. (22)], and Lieb 

[Eq. (24)]. The Sobolev inequality is given by Eq. (25). All 

quantities are dimensionless. 



(C/a) 


N 


Bargmann 
Limit 


Schwinger 
Limit 


Lieb 
Limit 


Sobolev 
Inequality 


1.65 


1 


2.19 


1.36 


1.49 


1.00 


1.68 


1 


2.25 


1.41 


1.53 




1.90 


1 


2.76 


1.81 


1.84 




1.98 


1 


2.95 


1.96 


1.96 




2.00 


1 


3.00 


2.00 


1.99 




6.00 


1 


21.00 


18.00 


10.33 




6.45 


2 


24.03 


20.80 


11.52 




7.00 


2 


28.00 


24.50 


13.02 




8.50 


2 


40.38 


36.13 


17.43 




8.92 


5 


44.24 


39.78 


18.73 





Reference 5. 



' The relative strengths of the Lieb, Schwinger, and Bargmann limits are depend- 
ent upon the potential V(r). For example, the table analogous to table 1 for the attrac- 
tive three-dimensional square well shows that the Bargmann limit is best in the range 
2.47 < Foo* < 10.74 and the Lieb limit is best when V<fi? > 10.74, where F is the 



depth of the well and a is the width of the well. The Sobolev limit states that bound 
states do not occur when V<fl < 2.11 whereas the exact numerical result states that 
bound states do not occur when V^f < 2.47. 
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8. Appendix 

In this appendix, we outline the evaluation of the integral that appears in eq. (21); namely, 

hi = Jo (1 +VT 2 In [1 +J?)/(1 -v)\ dri . (Al) 

We let 1 + ?/ = w and write 

2 

hx = \t w ' 2 \Hw) ~ In (2 - w)} dw . (A2) 

Referring to integral tables, we find that 

T — 1 _i_ Lim (2— w) i /o \ 1 ,.„x 

hl ~ Y + ,-2 V /n(2 " w) = T* (A3) 
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